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We demonstrate an unexpectedly strong surface-plasmonic absorption at the interface of silver 
and liigh-index dielectrics based on electron and photon spectroscopy. The measured bandwidth 
and intensity of absorption deviate significantly from the classical theory. Our density-functional 
calculation well predicts the occurrence of this phenomenon. It reveals that due to the low metal- 
to-dielectric work function at such interfaces, conduction electrons can display a drastic quantum 
spillover, causing the interfacial electron-hole pair production to dominate the decay of surface 
plasmons. This finding can be of fundamental importance in understanding and designing quantum 
nano-plasmonic devices that utilize noble metals and high-index dielectrics. 

PACS numbers: 42.79.Wc, 73.20.Mf, 78.20.-e, 78.68.+m 


Surface plasmons (SPs), collective oscillations of 
conduction electrons at a metal-dielectric interface, 
have attracted interest for several decades [l|-Q ■ 
Nanomaterials that strongly absorb visible light through 
plasmonic effects could be very important for solar- 
energy devices 0H. It is normally assumed that 
classical theory, with prescribed frequency-dependent 
bulk permittivities, reliably captures the SP properties, 
while quantum effects, despite their academic interest @- 


14], usually have negligible effect in practical systems. 


Here we show that contrary to conventional wisdom, 
quantum effects can play a crucial role for SPs at 
the interface of silver (Ag) and (practically any) high- 
index dielectrics. Our density-functional calculation 
and spectroscopic measurement indicate the existence of 
remarkable non-classical plasmonic absorption in such 
systems. Interfacial electron-hole (e-h) pair production 
B Mi can become the predominant decay channel 
of SPs, exceeding the ordinary phonon scattering (the 
Drude loss) or interband transitions (the dielectric loss). 
The quantum origin of such plasmonic absorption has 
been largely overlooked in nano-plasmonics research. 

High-index dielectrics have been long considered as 
superior gate insulators in nano-electronics [20]. In the 
recent years, they have also attracted much interest in 
nanophotonics on the conversion of SPs to hot electrons 
B Sr 241- However, very limited attention has been paid 
to their quantum electronic properties that can reversely 
modify the SP response. The uniqueness of high-index 
dielectrics, compared with low-index dielectrics lies in 
their large electron affinity, high static permittivity, and 
thereby a much lowered work function to the conduction 
electrons in metal. This allows the electrons to undergo 
a deep quantum spillover into the high-index dielectrics 
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FIG. 1: Schematics for the coupling of surface plasmons and 
interfacial electron-hole pairs at the interface of a metal and 
a dielectric. The e-h pair production can be viewed as dipole 
transitions, across the low barrier of the high-index dielectric, 
driven by the out-of-plane electric field of surface plasmons. 


extending beyond the Thomas-Fermi screening length. 
The in-plane SP propagation can strongly couple with 
the out-of-plane e-h pair production, as illustrated in 
Fig. [T| This process is insignificant at a metal-low-index 
interface [2, HR 16 1 due to the high barrier and thus 
insufficient quantum spillover. But it is remarkable at 
a metal-high-index interface and can lead to intensified 
energy dissipation and widened absorption spectrum. 


Our theoretical study utilizes a generalized jellium 
density-functional model that incorporates the crucial 
properties of Ag and dielectrics. A uniform Ag slab with 
a thickness d = 100 ao (ao is the Bohr radius) along z 
is clamped by a dielectric on both sides (see the inset of 
Fig. [2]). The conduction electrons are governed by the 
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Kohn-Sham and generalized Poisson equations [25 
2 to dz 2 


+ V eii [n_(z)\ <p v (z) = e v ip v (z), (1) 

_d_ 

'dz 


— \ e{z\uj)\ UJ ^ 0 —^{z) =4tt e{n_(z) -n + (z)}. (2) 


ipv(z) and e v are the eigenfunctions and eigenvalues. 
n_{z) is the electron number density and is a sum of 
| </?!,(£) | 2 over the occupied orbitals up to the Fermi level 
/i F [281, l29(. n + (z) is the space-dependent positive-jellium 
density; e(z; u>) is the space- and frequency-dependent 
background permittivity that accounts for the screening 
from valance electrons. The effective potential reads 

V eH [n_{z)\ = -e<P{z) + V xc [n_(z)\ + a(z), (3) 



Position z (A) 


in which V xc [n_(z)\ is the exchange-correlation potential, 
and a(z) is the space-dependent electron affinity in a 
flat-band picture [30|. A generalized Green’s function 
Q(z,z',q; cj) for this system is deduced in conjunction 
with the Poisson equation for the statics (the in-plane 
wavenumber q 0 and frequency w —> 0 ) and dynamic 
response. (See Supplemental Material for details.) 

Several high-index dielectrics, AI 2 O 3 , HfC> 2 , and Ti0 2 , 
are investigated; and the common medium-/low-index 
dielectrics, Si0 2 and air, are included for comparison. 
Table [J lists the material properties adopted or obtained 
in our calculation. Figure [2] shows the near-interface 
ground-state electron density profile n{z)/n and effective 
potential I4ff(2)- With increasing static permittivity and 
electron affinity, the potential barrier (work function) 
drops by as much as 2 eV, and an increasing number 
of electrons spill from Ag into the dielectrics. This 
behavior can be quantified by a characteristic spillover 
depth, £ = J 0 -oc dz n_(z)/n representing the distance 
up to which the spilled density would extend if it had 
the constant bulk value n [2S, [l 6 j. As shown in Table Q] £ 
of HfC >2 or TiC >2 is 2 to 3 times greater than that of air, 
and has approached the Thomas-Fermi screening length, 
l T f « 0.58 A of Ag. The actual density tail penetrates 
several times further, as displayed in Fig. [2j 


TABLE I: Adopted material constants of the dielectrics 
[I3,[32-H: electronic bandgap E g , electron affinity a D , static 
permittivity e D o = £d(<^ 0), and permittivity at 532 nm 

optical frequency £ 0 (^ 532 ). Calculated Ag-to-dielectric work 
function W and spillover depth £. 


Material 

E g (eV) 

Qd (eV) 

Cdo 

Cd (^532) 

W (eV) C (A) 

Air 

“00” 

0 

1 

1 

4.26 

0.18 

Si0 2 

9.0 

0.95 

3.9 

2.13 

3.22 

0.24 

AI 2 O 3 

OO 

bo 

1.70 

9.0 

3.13 

2.73 

0.32 

Hf0 2 

5.8 

2.14 

25.0 

4.33 

2.47 

0.39 

Ti0 2 

3.2 

3.00 

86.0 

6.55 

2.26 

0.58 


FIG. 2: Calculated ground-state electron-density profiles, and 
effective-potential profiles with respect to the Fermi level /r F • 


To see how quantum spillover influences the dynamic 
excitations, we implement a time-dependent linear- 
response calculation [13, SI]. The system-mediated 
effective interaction between two external charge sheets 
at Z and Z' is 371 


W(Z, Z',q-,u) = J dzdz' Q(Z, z, q;u) 


(4) 


x x(z > z',q;u)g(z > ,Z',q;u), 


where x(z, z', q; u> ) is the susceptibility. (See Supplemen¬ 
tal Material.) The surface response function, g(q;ui) = 
(g/27r)e D (o;)Im>V(^0j Z 0 , q ; w) describes the amplitude of 
surface excitations caused by an external charge sheet at 
Zq, which we take to be 50 ao outside the Ag 0. Fig¬ 
ure [31(a) gives the calculated g(q;u>) at a representative 
wavenumber q = 0.05 aj/ 1 , which amounts to a 6.6 nm 
wavelength, a typical length scale relevant to both optical 
and electronic excitations. The main peaks in Fig. [31(a) 
correspond to the SP excitations. As the dielectric index 
rises, they change from a narrow 3.57 eV resonance for 
air to a broad feature around 2.51 eV for TiC> 2 . Fig. [31(b) 
shows the induced density variation corresponding to the 
peak frequencies, and confirms the surface-mode profiles. 

The spectral broadening (and the associated SP 
damping) here takes place without phonon scattering or 
interband transitions (our density-functional model does 
not include these mechanisms). The intrinsic dissipation 
occurs via interfacial e-h pair production. According to 
Feibelman et al. [3, ll 6 ], the decay rate of SPs through 
e-h pairs is generally proportional to the imaginary part 
of a d ± parameter, which is subsequently proportional 
to the spillover depth ( through sum rules. A classical 
system with £ = 0 cannot produce any out-of-plane e- 
h pairs. For metal-low-index contacts, ( is practically 
too small to activate this dissipation channel. Only for 
the metal-high-index contacts studied here, ( becomes 
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FIG. 3: (a) Calculated dynamic surface response function, 
(b) Calculated induced density variation (scaled by e D for 
clarity) of surface plasmons at peak response energy, (c) 
Calculated energy-momentum loss spectrum r(q,w) for Ag- 
Hf02, plotted in logarithmic scale. Red means high loss and 
blue means low loss. Note that for a 100 aq thin Ag slab in 
our model, there can be more than one bulk-plasmon curves. 


appreciable and thus e-h pair production can induce 
significant effects. 

To see the underlying physics comprehensively, we 
calculate the energy-momentum loss spectrum / 1 (q. uj ) 
for a high-energy electron penetrating the system 37l.l38|, 


r(q,uj) = —C J dZdZ' cos 
x lmW(Z, Z',q;uj 


-(Z-Z’) 

Via 


(5) 


where C is a universal constant and Uj n is the incoming 
velocity of the electron (80 keV in our calculation and 
experiment). For a conventional Ag-vacuum contact 
(the work function W = 4.26 eV), the spectrum should 
contain a clear Lindhard boundary between the collective 
excitations (surface and bulk plamsons) and the bulk 
e-h pair production jssj]. Unless the Landau damping 
comes in at extremely large momenta, SPs cannot decay 
into bulk e-h pairs without violating energy-momentum 
conservation (see Supplemental Material for details). 
However, our calculated spectrum Fig. Etc) for an Ag- 
HfCU contact {W = 2.47 eV) exhibits a thoroughly 
blurred region between the collective excitations and the 
bulk e-h pair production. This is the region of interfacial 
e-h pair production. For a metal-dielectric contact of a 


small work function (low barrier), electrons can actively 
move out-of-plane irrespective of how small the in-plane 
q is. SPs can spontaneously decay into interfacial e-h 
pairs without violating energy-momentum conservation. 
This manifests the importance of ground-state properties 
to the optical-frequency behaviors in such systems. 

Looking for experimental evidence, we perform 
electron energy-loss spectroscopy (EELS) measurement 
(aberration corrected Zeiss Libra transmission electron 
microscope (TEM)). The acceleration voltage of the 
microscope is 80 kV, which permits an excellent energy 
resolution of 120 meV (FWHM of the zero-loss peak). We 
e-beam evaporate 20 nrn Ag and 20 nm varied dielectrics 
sequentially onto carbon-supported TEM grids. Figure 
[4] shows the obtained energy-loss spectra. The broad 
peaks in the range of 1.5 to 1.8 eV for all dielectrics come 
from_the contact between 5 nm amorphous carbon and 
Ag 0 - These peaks are not of interest. The sharp 
peaks around 3.8 eV for all dielectrics are the bulk- 
plasmon resonances inside Ag 41]. The varied peaks 
consistently moving from about 3.5 eV for SiCU to about 
2.8 eV for TiC >2 are the SP resonances at the Ag-dielectric 
interfaces. There is a clear trend that, with increasing 
index of the dielectrics, the SP peak-width gradually 
broadens. This trend signifies an enhanced damping 
and shortened lifetime of the SPs, in agreement with 
our theoretical prediction in Fig. [3] By comparison, 
the bulk-plasmon resonances, which decay primarily 
through phonon scattering (the Drude loss) and bulk e-h 
pair production (only at extremely short wavelengths), 
are much less affected by the varied dielectrics. This 
observation testifies to the dominant role of interfacial 
e-h pair production, as opposed to phonon scattering, in 
the SP absorption at Ag-high-index interfaces. 
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FIG. 4: Measured electron energy-loss spectra. The colored 
column plots display the electron counts. The colored curves 
show the peaks via the standard multi-Lorentzian fitting 
(after subtracting the zero-loss peak). 


In order to directly connect to optical measurement, 
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we further develop a mesoscopic formulation to quantify 
the crossover from the phonon-dominant regime to the e- 
h-pair-dominant regime. The p- wave reflection coefficient 
at a metal-dielectric interface allowing quantum spillover 
can be derived by generalizing the original formulation 
by Feibelman et al. [2, ll5Hl^ , 


r {p \q,u) 


€p 

e D 


hvM 
K M 


e M 




6 m ^D 


[' d±q 2 + d||K D K M ] 

-■ ( 6 ) 

[d ± q 2 - d n K D K M \ 


Here, e M (w) = e B (to) — w 2 /(w 2 + i 7 P w) is the Lorentz- 
Drude permittivity of bulk Ag (see Supplemental 
Material), /c D (g,w) = \Jq 1 — e D (w)u; 2 /c 2 , re M (g,w) = 
\Jq 2 — e M (w)w 2 /c 2 are evanescent wavenumbers. cC and 
d|| are two complex-valued parameters, associated with 
the long-wavelength electron-density oscillation near the 
interface. If oC and dy both vanish, then Eq. (©) 
is reduced to the standard p-wave Fresnel reflection 
coefficient that embodies classical SP resonance when the 
denominator ^ is near zero, dy typically vanishes 
for uncharged surfaces. The real part of d ± gives the 
nonlocal correction to the SP dispersion. The imaginary 
part gives the more important e-h pair loss strength. Our 
density-functional calculation gives Imdj_ ~ 1.1 A for 
Si 02 , 1.5 A for AI 2 O 3 , 2.1 A for Hf 02 and 3.0 A for 
Ti 02 0. Note that Imd x is larger than the ground-state 
spillover depth ( listed in Table HJ 

The phonon (Drude) loss of Ag is small compared 
with visible frequencies, so a linear expansion to j p /co in 
Eq. © around the conventional SP dispersion curve ui(q) 
(solved from q 2 = [e M e D /(e M + e D )]w 2 /c 2 ) is legitimated. 
Here e M (w) is the q p —> 0 limit of e M (w). We can 
determine a dissipation ratio rj[ui(q)] (see Supplemental 
Material) that evaluates the strength of e-h pair loss 
versus the phonon loss at each SP-resonance frequency 
along the dispersion curve ui{q), 


V[u(q)\ 


2 |e M (w )| 2 w 4 Imd_L 
VIcmMI -e D (w) 7 P c ' 


(7) 


In Fig. ©a), we plot 77 as a function of ui for 
all the dielectrics, with the value of q p fitted from 
real Ag, and with Imd x varying from 1 A to 3 A. 
The red dashed-dotted line separates the phonon- 
dominant regime (r)(to) < 1 ) from the e-h-pair-dominant 
regime (??(w) > 1). The conventional SP-dispersion 
asymptotically approaches the frequency where |e M (w) | = 
e D (w) (the black dashed lines in Fig. ©a)). At this 
frequency (corresponding to very large q), the e-h pair 
loss always dominates, even for a medium-index dielectric 
like SiC> 2 . However, towards smaller frequencies (longer 
2-kc/u but still fairly large q), Ag-SiC >2 with Imd x ss 
1.1 A immediately drops to the Drude-dominant regime. 
Ag-Al 2 C >3 with Imdj^ « 1.5 A is at the onset where e- 
h pair loss starts to take effect in a wider range up to 


Surface-Plasmon Resonance Energy hui (eV) 
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FIG. 5: (a) Calculated dissipation ratio of the e-h-pair loss 
versus the phonon loss around the conventional SP dispersion 
curve, (b) Measured ultraviolet-visible reflection spectra, and 
calculated reflection spectra using classical roughness theory. 
The shaded region (A < 350 nm) is where silver approaches 
bulk-plasmon absorption and can have larger uncertainties. 


2 - 7 T c/uj ~ 400 nm or so. Ag-HfC >2 with Imd^ ss 2.1 A has 
an even wider e-h-pair-dominant range beyond 500 nm. 
Ag-TiC >2 with Imd^ ~ 3.0 A is fully e-h-pair-dominant 
covering the entire visible spectrum. This result suggests 
possible observation, via optical approach, of strong 
and broadened resonant absorption around the large- 
momentum SP frequency (around the black dotted-lines 
in Fig. ©a)), especially for Ag-high-index contacts. 

We perform ultraviolet-visible (UV-Vis) spectropho¬ 
tometry measurement (Varian Cary 500). We e-beam 
evaporate 20 nm Ag films onto ultra-smooth (rough¬ 
ness 6 < 3 A) z-cut single-crystal quartz substrates, 
and conformally coat 20 nm varied dielectrics on the 
top via atomic layer deposition (ALD). The 20 nm Ag 
on the dielectric side has a roughness S ~ 2 nm and 
a correlation length a ps 35 nm measured by atomic- 
force microscopy (AFM). This slightly roughened sur¬ 
face provides the desired large momenta to excite SPs 
under normally incident light without a need of prism 
coupling 0, rilffit For convenience, we deposit addi¬ 
tional 200 nm Ag on the bottom of the substrates to 
block the transmission, so we can measure the reflection 
spectrum from the top (refer to the inset of Fig. ©b)). 
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In the wavelength range of 350 to 700 nrn, the quartz 
substrates are completely transparent and the 200 nm 
Ag on the bottom serves as an ideal mirror. Any drop 
on the reflection spectrum is due to the absorption from 
the 20 nm dielectric and 20 nm Ag on the top (diffu¬ 
sive scattering is negligible at this roughness). The solid 
curves in Fig- Etb) show the observed resonant absorption 
around the large-momentum SP-frequency of each dielec¬ 
tric. The reflection dip is not prominent for SiC> 2 , but 
continuously intensifies and widens from AI 2 O 3 to Ti 02 , 
as anticipated from Fig. (51(a). The observed absorption 
is not caused by interband transitions (dielectric loss) 
inside the dielectrics, as we have done separate transmis¬ 
sion measurement for the dielectrics on quartz alone (no 
Ag) to confirm their lossless nature in the wavelength 
range of interest (see Supplemental Material). It is not 
solely due to the phonon-scattering (Drude loss) in Ag 
either, as we have carried out a classical transfer-matrix 
calculation combined with the statistical-roughness the¬ 
ory following Kretschmann et al. jiil, E3 |. This calcu¬ 
lation encloses all the Drude absorption at the interface 
and bulk of Ag. The dashed curves in Fig. [5][b) show 
the reflection dips from the classical theory. They are 
sharper and located at longer wavelengths (determined 
by the er ss 35 nm correlation length), in contrast to 
the experimental curves. Only Si 02 exhibits a similar¬ 
ity between the experimental and (classical) theoretical 
curves. All the high-index dielectrics show large discrep¬ 
ancies. Based on our systematic theoretical development 
above, the e-h-pair-dominant loss is the most probable 
reason for the observed strong and broadened absorption 
of short-wavelength SPs at Ag-high-index interfaces. 

To summarize, we find that high-index dielectrics 
contacting with silver can exhibit enhanced surface- 
plasmon absorption due to the quantum-spillover 
supported interfacial electron-hole pair production. 
The quantum-electronic control of the static dielectric 
environment to the optical excitations on metal surfaces 
can bring on new applications in nanoscale light 
confinement and new insights in surface-plasmon to hot- 
electron conversion. 
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